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First-principles generation of Stereographic Maps 
for high-field magnetoresistance in normal metals: 
an application to Au and Ag. 
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Cittadella Universitaria, 1-09042 Monserrato (CA), Italy 

(Dated: February 2, 2008) 

About thirty high-field magnetoresistance Stereographic Maps have been measured for metals 
between Fifties and Seventies but no way was known till now to compare these complex experimental 
data with first-principles computations. We present here the method we developed to generate 
Stereographic Maps directly from a metal's Fermi Surface, based on the Lifshitz model and the 
recent advances by S.P. Novikov and his pupils. As an application, we test the method with an 
interesting toy model and then with Au and Ag. 

PACS numbers: 72.15.Gd,02.40.Re,03.65.Sq 



I. INTRODUCTION 

The relevance of the geometry and topology of the 
Fermi Surface (FS) in physical phenomena is well known 
since Thirties, when Justi and Scheffers showed evidences 
that the Fermi Surface of Gold is operi. 

One of the most striking examples of such phenomena 
is the behaviour of magnetoresistance in monocrystals at 
low temperatures in high magnetic fields. Evidences for 
this effect was first discovered theoretically by Lifshitz 
and his Karkov school in Fifties^ by studying galvano- 
metric effects in metals without any special assumption 
for the electron energy spectrum, and it was verified ex- 
perimentally by Gaidukov et al. shortly afterwards^'^. 

What was clear from those early works is that, under 
the conditions stated above (and apart from the excep- 
tional case when the density of electrons and holes are 
equal) the magnetoresistance behaviour is dictated only 
by the topological properties of orbits of quasi-momenta 
(see fig^): as the magnetic field H grows, the magnetore- 
sistance saturates isotropically to an asymptotic value if 
the orbits are all closed, while it grows quadratically with 
H if there are open orbits; moreover, in this last case the 
magnetoresistance is not isotropic and the conducibility 
tensor a has rank 1. 

Many works about this effect were pub- 
lished between Fifties and Seventies from 
both the experimenta l^i^i^i^i'^1^'^'^^'^^'^^ and 
theoretical^i^'^i^^i^-^^^i^^i^^i^^i^"i^^ point of view; in 
particular, Stereographic Maps (SM) were experimen- 
tally built by plotting in a stereographic projection all 
magnetic field directions H in which a quadratic rise of 
a was observed (see figEJ. In those times the interest 
on these magnetoresistance effects was due mainly to 
their utility as a tool to determine FS properties rather 
than as phenomena in their own right, and in particular 
SM maps provided information about FS topology: e.g. 
clearly if a grows quadratically for some direction of H 
then the FS must be open, and further analysis can lead 



to discover the directions of the openings. 

Between Fifties and Seventies SM were experimentally 
found for about thirty metals but, despite the theoret- 
ical efforts, no way was found to generate them with 
first-principles calculations and therefore no accurate di- 
rect verification of the Lifshitz model is available to 
date, except for the qualitative sketches by Lifshitz and 
Peschanskiii^ (see figEl; in particular it was not known 
till now how closely the semiclassical model is able to 
reproduce these complex experimental data and whether 
further purely quanto-mechanical corrections are needed. 

As the magnetoresistance methods were replaced by 
newer and more accurate tools to study FS, no new SM 
was experimentally produced and the problem was even- 
tually abandoned; it was only in Nineties that the beau- 
tiful topological structure underlying this phenomenon, 
this time considered in its own right rather than as a 
tool for something else, was fully discovered, making this 
way finally possible the construction of an algorithm able 
to reproduce from first-principles the experimental data 
about the dependence of a on the direction of H for any 
given Fermi Function (FF) £. 

In 1982 indeed S.P. Novikov^^ recognized the purely 
topological character of the problem and later^^ ex- 
tracted the following generic picture from the work of 
his students A. Zorich^^ and I.V. Dynnikov^^'^^'^^: once 
a Fermi Surface (FS) is given, if open orbits arise for 
electrons quasi-momenta for some direction of the mag- 
netic field, then the set of such directions are sorted in 
some finite number of "islands" (e.g see fig. [2)); to each 
of these islands it is associated a new quantum invariant 
L (an irreducible Miller index of the lattice) that defines 
the dynamic of the semiclassical system in the following 
way: each open orbit corresponding to a ii" that belongs 
to an island labeled by L is a finite deformation of a 
straight line parallel to the vector product H x L. All 
H that do not fall in any of these islands give rise only 
to closed orbits, except for a negligible set of exceptional 
directions that we will disregard here. 

In this article we present the method, suggested to us 
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FIG. 1: Behaviour oi p — cr~^ \\i metals with closed and open FS^^. (Closed) p is isotropic and saturates unless the density of 
electrons and holes coincide, in which case p ~ . (Open) p is highly anisotropic and it shows qualitatively different behaviour 
in minima and maxima (resp. Q2 and Q\ in the picture): in maxima p ^ , in minima it saturates {0 is the angle between H 
and the crystallographic axis). 



by L Dynnikov, that we implemented, based on this pic- 
ture, to detect the directions of the magnetic field for 
which open orbits of quasi- moment a appear. This algo- 
rithm allows us to predict the SM of a metal according 
to the Lifshitz model, namely to determine the part of 
the "islands" that the semiclassical approximation is suf- 
ficient to detect. As an application, we first study the SM 
of a very rich and simple toy model and then we gener- 
ate SM for Au and Ag and compare our results with the 
ones obtained experimentally more than forty years ago 
by Gaidukov (and never repeated since then); this repre- 
sents the first rigorous check of Gaidukov's results from 
first-principles and shows that, even though the SM were 
measured just at the lower threshold for the semiclassi- 
cal approximation to hold, the Lifshitz model is able to 
reproduce rather accurately the experimental data. 



II. OVERVIEW OF THE LIFSHITZ MODEL 

A great deal of energies have been invested in Fifties 
ans Sixties in the study of magnetoresistance in a metal. 
After the discovery by Kapitza^^ of a linear increase of 
the resistivity with a magnetic field in a number of met- 
als, it had be shown by Justi^ that metals could be di- 
vided in two categories: for the first one the resistivity 
saturates with the increase of the magnetic field (e.g. in 
Cu, Na, Al) while for the second one it grows quadrati- 
cally (e.g. in noble metals). 

It was Peierls^^ the first to recognize that this anoma- 
lous behaviour was due to the departures from the free- 
electron model, but it was especially thanks to the 
theoretical works of I.M. Lifshitz and his schoolS&i^ 
and to the experimental results of Alekseevskii and 
Gaidukov^'^'^'^'^'^'^'^^'^^ that a systematic and thorough 
study of this phenomenon was carried on and fully re- 
vealed its critical dependence on the topology of the FS. 
In particular, Lifshitz was the first to study the system 
in its full generality, making no assumption on the form 
of the FF. 



Let us review quickly the mathematics of the model to 
show the role topology has in it^^'^^. Since no analytical 
method to solve exactly the Schrodinger equation under 
a generic periodic potential is known, we will make, as 
usual, a first approximation introducing the semiclassi- 
cal model, i.e. we neglect the electron-electron interac- 
tion and consider the motion of a single electron in an 
infinite crystal. According to the semiclassical approxi- 
mation, the Schrodinger equation for a single electron in 
a three-dimensional lattice F with a F-invariant poten- 
tial in presence of an electric field E and a magnetic field 
H reduces to the following (semi-) classical equations of 
motion: 



p=-e [E ^ ^^Vn{p) X H] 

where £n{p) is the energy function for the electron occu- 
pying the band n. 

Physical constraints limit the range of the magnetic 
field: electrons become aware of the FS topology only if 
the mean free path is long enough to traverse considerable 
portions of it, that requires ujc^ 1 (and so a field H > 
lOT), a pure crystal and very low temperatures; on the 
other side, to avoid magnetic breakdown we also must 
have H < lO^T. 

The semiclassical approximation makes the problem 
look at first sight as a standard classical mechanics sys- 
tem but it turns out that it is instead deeply different 
from all of them. The difference is not in the analytical 
expression of the equations, that is evidently the same, 
but rather in the topology of the phase space: in classi- 
cal mechanics indeed the momenta belongs always to a 
linear space, no matter how complicated the base space 
is; here instead the base space is topologically trivial (the 
whole three-space) but the momenta are triply periodic 
by the Bloch theorem, since we identify all momenta that 
differ by a vector of the dual lattice F*. In more rigorous 
terms, while in classical mechanics momenta would be- 
long to the linear space R^, in this case they are defined 
only modulo a vector of the reciprocal lattice and there- 
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fore they belong to the first Brniouin zone, i.e., in other 
words, to the three-torus = R^/F*. 

This difference is essential because it is exactly what 
brings topology (in this case "periodic topology" ) in play: 
the three-torus, differently from the three-space, has a 
non-trivial topology and its presence has a strong influ- 
ence on the dynamics of the system. Indeed, if both E 
and H are constant, the pair of (systems of) equations 
de-couples and the problem reduces to study the orbits 
of the quasi-momenta in the first Brillouin zone (with the 
obvious boundary conditions dictated by the periodicity) 
under the equation: 



p = —e 



' ldS{p) _ 
X H 

c op 



As often happens, topological effects are due to the 
magnetic field rather than to the electric one and there- 
fore we can safely put E = and e = c = 1 and we are 
finally left with the equation: 



d£{p) 

dp 



xH = {p,S{p)}, 



where {, }jj is the so-called "magnetic bracket": 

It is well known that this system is over-integrable since 
it has two integrals of motion, namely the Hamiltonian £ 
and the component p • H of the quasi-momentum in the 
magnetic field direction; nevertheless it was not fully un- 
derstood until Eighties the relevance of the fact that the 
latter integral is not a well defined singie-valued function 
in but rather a muiti-valued function, exactly in the 
same way the angle is just a multivalued function in 
the circle S^. 

Let us point out that not even the most elementary 
systems with multivalued first integrals have been ob- 
ject of study till recent years by the dynamical systems 
community, probably because no such system arises nat- 
urally from classical mechanics problems. This appear- 
ance of a multivalued first integral, whose occurrence is 
due to the non-trivial topology of and could not arise 
in a topologically trivial space like R^, it's enough to 
transform an otherwise trivial dynamical system in an 
extremely rich one. For example, cutting level surfaces 
of £ with level surfaces of another well-defined function 
on would lead only to closed orbits in the Brillouin 
zone, while in our case open orbits do generically arise, 
orbits that may in principle fill the whole FS in the same 
way a straight line with 3-irrational slope would fill the 
whole Brillouin zone. 

The topology of the orbits, namely whether they are 
open or closed (in the repeated zone scheme), is not just 
a mathematical curiosity: it has been first showed by Lif- 
shitz indeed that the cause of the quadratic growth of the 
magnetoresistance is exactly the presence of open orbits. 
A rigorous deduction of this fact can be made using the 
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FIG. 2: (a) SM of Gold measured by Gaidukov in 1959^ (b) 
Qualitative sketch of the previous picture obtained by Lifshitz 
from his topological analysis^ 



Boltzman transport equation^^ but simpler justifications 
of it can be made^^ at a phenomenological level based on 
the Einstein relation a Ve^n/£F^ where V is the diffu- 
sion coefficient and £f the Fermi Energy. When all orbits 
are closed indeed the diffusion in the plane perpendicular 
to H consists in jumps by an amount of the order of the 
cyclotron radius vh = cpp/eH with frequency ~ 1/r, so 
that V Iv ^h/^ ^ — <^o/(^ct)^; now suppose 
instead that quasi-momenta open orbits appear instead, 
say in the x direction: then the electrons will move along 
the y direction, since the orbits in the real space are ro- 
tated by 7r/2 with respect to orbits in the momentum 
space, so that in the x direction V has more or less the 
same H dependence found for the closed orbits but in 
the orthogonal direction electrons move as free particles, 
pT, and therefore in this case a is not 



yy 



namely V, 

isotropic anymore and axx — (^^/{'^c^Y^ (^yy — c^o- 

The theoretical results by Lifshitz and his school urged 
Alekseevskii and Gaidukov to start conducting careful 
experiments that turned out to be in perfect agreement 
with the model. The most interesting experimental result 
for us is the stereographic projection of the special direc- 
tion for the magnetic field (Stereographic Map), namely 
the map on the unitary disc that shows for which direc- 
tions the magnetoresistance grows quadratically. Such 
map indeed, in the semiclassical approximation, depends 
solely on the orbits topology and therefore it is totally 
determined once a Fermi Surface is given. A great ef- 
fort has been spent to find some kind of algorithm able 
to produce this map from a generic FSiSiiLii, but the 
topological tools to solve the problem were not known to 
physicists at that time and eventually the problem was 
left unsolved. 



III. OVERVIEW OF THE RECENT 
TOPOLOGICAL RESULTS 

The interest in this problem revived in 1982 when 
S.P. Novikov found out that it was a perfectly suited 
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(c) ^ ' (d) 

FIG. 3: FS with rank 3 (a), 2 (b), 1 (c) and (d). 



case where to apply his newly introduced Morse- Novikov 
theory^^ , namely the study of the topology of level sur- 
faces of multivalued functions. 

Fundamental results were found in Eighties and 
Nineties by his pupils A.V. Zorich and LA. Dynnikov and 
from them Novikov later extracted the following picture: 
once a "complicated enough" Fermi Function is given (we 
will clarify this concept below), a fractal is determined on 
the set of the stereographic projections of the magnetic 
field directions; the fractal consists of smooth polygons 
(islands, or "stability zones") that contain all possible 
lattice directions and generically meet each other in a fi- 
nite number of points; every such polygon is labeled by a 
Miller index L and moreover to every magnetic field di- 
rection are associated two values of the energy em,M{H). 

The meaning of this picture is the following. Let £ 
and Ef he the Fermi Function and Fermi Energy of a 
metal and suppose that we want to know the asymptotic 
behavior of trajectories of quasi- momenta for some mag- 
netic field H whose stereographic projection lies inside 
a polygon labeled by the Miller index L: the answer is 
that if Ep lies between em{H) an eM{H) then open or- 
bits exist and moreover they are all finite deformations 
of a straight line with direction H x L; if Ep is instead 
smaller than em or bigger than eu then all orbits are 
closed. 

The directions for which em{H) = eM{H) are non- 
generic and Novikov conjectured that they form a set of 
measure zero in the disc. The orbits of quasi-momenta 
induced by magnetic fields with those directions are much 
more complicated than the generic ones and for each such 
fixed direction there is only one energy value for which 
open orbits appear, while at every other energy all orbits 
are closed. Nevertheless this case may be relevant since 
it may explain deviations from the quadratic growth law 
detected in several metals for special directions of the 
magnetic field and it is currently under investigation by 
Novikov himself and Ya. Maltsev^^. 

Below we review in some detail the topological ideas 
on which the method is based. The first important fact 
is that^^ every FS is topologically equivalent to a sphere 



with some finite number of handles attached (such num- 
ber is called the genus of the FS). It is an intuitive fact 
that we need a genus bigger than zero to have an open FS, 
but this condition is not sufficient: it is enough to think 
to a sphere in the center of the first Brillouin zone and 
to add to it handles without ever touching the boundary 
of the zone. 

The concept that detects whether a FS is open or not 
is its topological rank, namely the complement to three 
of the biggest number of "linearly independent" pairs of 
planes (i.e. planes whose perpendiculars are linearly in- 
dependent) that can enclose the FS. Examples of FS of 
rank from three to zero are shown in figO for example 
a sphere can be enclosed between three linearly indepen- 
dent pairs of planes, while a cylinder can be enclosed only 
between at most two of such pairs and so on. The case of 
interest for the phenomenon in study is the rank-3 one, 
since in the other ones either all orbits are closed (rank 
zero) or generically closed (rank one) or open but with 
an obvious asymptotic direction (rank two). 

Despite its triviality, genus-2 rank-2 case turns out 
to be paradigmatic and it is worth describing it in de- 
tail. Examples of such surface are shown in figEllSlb)^! 
namely a pair of parallel rectangles joined by a cylinder 
(incidentally, this is exactly the topology of the FS of 
Tin). The normal to the rectangles is a lattice direction 
that we will call L. 

Suppose first that H induces closed orbits on the cylin- 
der (as in fig^: then H induces also open orbits, since 
any orbit that does not lie entirely on the cylinder will 
be bound to stay on only one of the two plane sheets 
and therefore will be open (it has no way to turn back). 
Moreover, the surface is clearly enclosed between a pair 
of parallel lattice planes and therefore the orbits lie en- 
tirely in a finite width strip and, as Dynnikov^^ showed, 
pass through it; the Miller index common to the two 
lattice planes is clearly L and it is the quantum invari- 
ant associated to H. Suppose now that there exist a 
plane perpendicular to H that intersects both bases of 
the cylinder, so that no closed orbit is cut on the cylin- 
der by any plane parallel to it: then, assuming H fully 
irrational to simplify the discussion, every orbit will be 
closed. Indeed take any point on the FS and follow its 




FIG. 4: A genus-2 rank-2 surface. For this choice of H closed 
loops are cut on the cylinder and therefore the critical saddles 
through pi and p2 are half-opened (figEIii)) and open orbits 
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FIG. 5: (a) Generic picture in case open orbits exist (here ^ = 8); (b) warped plane filled by open orbits (the two discs shown 
are bases for two cylinders of closed orbits^ (not shown in fig.) 



orbit in one of the directions: since H is fully irrational, 
at a certain point the plane generating the orbit will cut 
a copy of the cylinder (in the extended zone picture) on 
both bases, i.e. the orbit will turn back; since exactly the 
same happens in the other direction, the orbit is closed. 

Our analysis shows that the SM corresponding to such 
a FS, as it has been obtained experimentally for Tin^, 
contains a single island, whose boundary is given by the 
set of directions of planes that are tangent to both bases 
in diametrally opposite points. For example, if the topo- 
logical cylinder is actually a right circular cylinder of ra- 
dius r and height /i, the island is a circle with center in 
the L direction and radius h/r. The Miller index associ- 
ated to this island is of course the same L above and the 
direction of open orbits is given by the vector product 
LxH. 

Let us come now to the rank-3 case. A fundamental 
result by Zorich^^ is that, no matter how complicated 
the FS is, the one above is the generic behaviour for all 
directions "close enough" to rational, namely the open 
orbits, when they arise, lie on components of the FS that 
are exactly "warped planes" separated by cylinders of 
closed orbits. For each fixed direction, all such planes 
are parallel to each other and are in even number. The 
number of pairs is bound from above by ^/2, where g 
is the genus of the surface; the number of cylinders is 
bounded from below by — 1 (see figlSJ. The difference 
with the rank-2 case is that now the surface may be split 
in many different ways in "warped planes" and cylinders 
(even in infinitely many ways in special cases), so that 
more than one zone can appear (possibly each one with 
a different Miller index). 

The key-point here is the strong constraint determined 
by periodicity to the type of critical point that can be 
met in the sections of the FS by planes. Consider the 
simple case of H = (0,0,1): assuming that there are 
open orbits and that we are not in a degenerate case, 
all critical points will appear at different levels and they 
will be of one of the types shown in figlHl since we are 
interested in open orbits, case {Hi) is irrelevant for us 
and case {iv) is forbidden by the boundary onditions. Let 
us now follow one of the open orbits: if it never meets 



any other curve, then after a period it meets itself again, 
meaning that the whole FS is just a single warped plane 
and everything is trivial (it is a genus- 1 surface with rank 
1 or 2); if it meets another orbit, then either it will be 
another open orbit (case (i)) and the two will annihilate 
each other (or, from another point of view, the open orbit 
will hit a closed one and bounce back), or it will be a 
closed orbit (case (ii)) and so the open one will simply 
engulf it and go on till it will meet again itself. There 
are only two possible outcomes for this process: either 
a warped cylinder (genus- 1, rank-1) or a warped plane 
(genus-1, rank-2), in each case with some finite number 
of plane holes that are the basis of the cylinders of closed 
orbits that separate them from each other (see fig|S}. 

The reasonment above holds for any H pointing to 
a lattice direction, since we can bring every such H to 
(0,0, 1) with a coordinate change that leaves the lattice 
invariant, but it is impossible to extend it to the case 
of "irrational" H (i.e. not directed along a lattice di- 
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FIG. 6: Possible kind of open saddles for rational H: (i) 
Fully open; (ii) Half-open; (iii) Fully closed; (iv) Fully open 
(impossible because of the boundary conditions). Only (ii) 
plays a role in the generic case since (i) arises only for rational 
directions of H. 
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FIG. 7: Possible choice of canonical loops in a genus-2 sur- 
face. Any other non-trivial loop (i.e. a loop that cannot be 
continously shrunk to a point) drawn on the surface cuts at 
least one of the base loops, i.e. it is "linearly dependent" by 
them. 

rections), since in that case the orbits are not periodic 
anymore. Nevertheless, continuity helps us there: indeed 
the heights of the cylinders of closed orbits depend (at 
least) continuously on the direction H and therefore the 
cylinders will survive small perturbations and open or- 
bits will still appear and will be bound to lie on small 
deformations of the same warped planes. 

The Zorich observation reveals the existence of a quan- 
tum invariant of the system that had previously been 
missed and, at the same time, explains why magnetic 
field directions in SM are sorted in islands: small defor- 
mations of warped planes preserve open orbits and lead 
to the same Miller index L, so there is a whole island 
around each rational H giving rise to open orbits and 
the same L is associated to all directions inside the is- 
land. From the construction it is clear that the open 
orbits have direction given hy H x so that experimen- 
tally it is enough to determine the direction of open orbits 
for two directions of the island to determine uniquely L. 

This discovery though does not solve completely the 
problem, since the set of directions "close enough" to 
all rational directions may not, in principle, cover the 
whole the disc. In Nineties I. Dynnikov showed^^ that 
the Zorich picture does represent the generic behaviour of 
the system in the following sense: given a Fermi Function 
£{p) and a direction open orbits arise only for a closed 
interval of energies [em{H)^ eM{H)] and the structure of 
open orbits is the one described above when em{H) = 
cm^H)] in particular this means that, for every possible 
direction of there is at least one value e of the energy 
such that H induces open orbits on f (p) = e. 

Dynnikov's results also lead to a conclusion rather in- 
teresting and unexpected: SM taken at different energies 
are compatible (i.e. no more than one L is associated 
to every H) and when we plot them all at the same 
time, building some kind of "all energies" SM, there are 
only two possible cases: either just one zone is present^, 
and that zone fills the whole disc, or infinitely many 
zones appear, distributed in a fractal-like way (see figEJ. 
Even though we cannot freely change the Fermi Energy 
of a metal, these fractals may in principle play a role in 
Physics because, for some special class of functions, there 
is an energy level at which the standard SM coincides ex- 
actly with the "global SM" (an example is provided in 



sec fV)) . 

IV. THE METHOD 

The topological picture discovered by Zorich makes rel- 
atively simple to build an algorithm able to reconstruct 
the magnetoresistance map from any given FS. 

Indeed, from what we said in the section above, it is 
clear that to investigate the topology of the orbits there is 
no need to study their asymptotics but rather it is enough 
to evaluate the Miller index associated to every magnetic 
field direction. Moreover, because of the Zorich result, it 
is enough to study what happens in the case of rational 
directions; this changes qualitatively the nature of the 
numerical problem we have to deal with, since the FS 
sections by a Miller plane are periodic and therefore, in 
principle, the whole orbits can be numerically evaluated 
with the desired precision. 

To evaluate the Miller index we need to use elemen- 
tary homology properties of loops on surfaces. It is a 
well known fact in topology that the algebraic sum of the 
intersections between two loops {intersection number) on 
closed compact surfaces is a homology invariant, i.e. it 
does not change if a loop is deformed with continuity or 
if it is replaced with one that forms with it the boundary 
of some surface. Moreover, if we agree to consider as the 
same loop all loops that are homologous to each other, 
we can give a Z-linear structure to this set that turns out 
to be isomorphic to Z^^, i.e. it is (freely) generated by 
2g (classes of) loops (see figQ). In this setting, the inter- 
section number becomes a bilinear antisymmetric non- 
degenerate form (,), namely a symplectic structure on 
Z^^, and in any canonical base {o^i^ Pj)ij=i^,,^g we have 
that 

< ai,/3j >= 1 , < ai.aj >=< /3i,/3j >= . 

The key observation is that the loops that lie on the 
"warped planes" of open orbits have all intersection num- 
ber zero with the closed loops living on the cylinders, 
since in that case there is no intersection at all. To eval- 
uate the desired Miller index therefore it is enough to 
find out the homology classes of all loops that have zero 
intersection number with the loops that form the closed 
cylinder and finally find out the lattice direction that 
these loops have in the reciprocal lattice F"*": since all 
of them lie on "warped planes" homologous to the same 
lattice plane, the result will be a two-dimensional sublat- 
tice of F"^ that will automatically give the Miller index 
we looked for. 

We implemented this method in the following way: 
once a FS M and a direction H are given, first of all 
we identify somehow a base (better if canonical) for the 
homology loops on M and provide methods to evaluate 
the intersection number with respect to them; then we 
determine all saddle points of type {i) and {ii) that H 
induces on M and evaluate the intersection number of 
the closed loops of each of these saddles with all loops of 



FIG. 8: (a) Rough sketch of the stabihty zones for a genus— 3 
surface^; (b) Example of section of a genus— 3 surface; (c) 
Component of the "prison bars" surface; (d) Basic component 
of the "prison bars" surface with basic cycles shown on it (each 
repeated twice). 



FIG. 9: (a) SM for E = cos(27rp:^) + cos(27rp^) + cos(27rp^); 
(b) detail of the fractal in coordinates {Hx,Hy, 1); (c) bound- 
aries and Miller indexes of a few stability zones evaluated an- 
alytically; (d) comparison between numerical and analytical 
calculations. 



the base above. Once all these data have been collected 
(and after many compatibility checks have been made) 
we identify the space of non-trivial loops that do not in- 
tersect them, namely we find a base for the space of all 
loops of M that lie on the warped planes^^-. Finally, we 
evaluate the directions of those loops in the reciprocal 
lattice F*, that will give us a set of directions identifying 
a 2-dimensional sublattice of F* that in turn will give us 
the desired Miller index. 

Our software core is a C++ library called NTC 
(Novikov Torus Conjecture) built over the graphics li- 
brary VTKi2& that provides the basic calls to build iso- 
surfaces meshes and perform elementary geometric and 
topological operations. What our library adds to VTK 
is mainly the capability of dealing with "periodic geome- 
try" and critical slices (e.g. planes that cut a surfaces in a 
saddle point) and the capability of evaluating topological 
quantities like homology classes of loops. 

The library has been released under the GPL license 
and it is downladable at the address http: / / ntc.sf.net / , 



V. A TOY MODEL 

It is well known^^ that a first rough expression for topo- 
logically non-trivial FS can be quickly obtained through 
the tight-binding approximation. 

For a simple cubic crystal this approximation leads in 



the lowest order to the following FF?i: 

E{p) = cos(27rpa^) + cos(27rp2y) + cos(27rj92) • 

It turns out that this FF is the simplest non-trivial 
case for our system: indeed its level surfaces are either 
spheres, when E < —1 or E > 1, or genus-3 rank-3 
surfaces shaped as a a sort of "three-dimensional prison 
bars" (fig. El). 

Several efforts have been made in Sixties to understand 
the topology of orbits in this elementary case, but the 
production of SM did not get more detail than the rough 
sketch shown in figlHI^a). 

This function is a perfect "toy-model" for testing our 
algorithm since its analytical expression is so simple to 
make possible verifying many things analytically. More- 
over, it satisfies the following property: its level surfaces 
{S = c} and {S = —c} differ only by a translation, so 
that a magnetic field H that gives rise to open orbits at 
energy c does the same at the opposite energy, i.e. all 
energy intervals giving rise to open orbits have the form 
[—e{H)^e{H)]; open orbits therefore arise for every di- 
rection of H at the level e = 0, namely the SM at the zero 
level is identical to the "global SM". It is easy to verify 
(even by hand!) that this SM has more than a zone and 
therefore fsec iIII|) it has infinitely many zones distributed 
in a fractal-like way. Finally, at the zero-energy level it 
is possible to get a simple analytical expression for the 
saddle points as function of the magnetic field direction, 
an information that is enough to obtain, in principle, the 
analytical expression of the boundary for any island. 
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FIG. 10: Miller indices associated to the biggest zones of the 
SM in the unitary square (figlHl)- 



We produced the SM in the square < Hx^y < 1, 
Hz = 1 by evaluating the Miller index associated to ev- 
ery direction in a 10^ x 10^ equally spaced square lattice 
(figin^b)) and then obtained the whole SM by symmetry 
(figElJa)). In figEtc) are shown the boundaries of the 
biggest zones obtained directly from the analytical ex- 
pression of the saddle points; as it is possible to verify 
from figE^d), there is a perfect agreement between nu- 
merical and analytical data. Similar calculations made 
for a piecewise smooth quadratic function with the same 
symmetries of £ showed a similar agreement even at en- 
ergies different from zero. 

Finally, in figE3we list the Miller indices correspond- 
ing to the biggest 26 zones together with an estimate of 
their area in the unitary square. A thorough description 
of the study of this FS can be found in [DLOa]^^. 



VI. NUMERICAL EXPLORATION FOR GOLD 
AND SILVER 

It is a well-established result that Gold and Silver have 
a genus-4 FS: both FS are spheres with four handles 
directed along the four diagonals of the cube (e.g. see 
figda)). 

Extremely precise estimates of the FF values can now 
be achieved through numerical calculations, for exam- 
ple using advanced tight-binding methods^^, but sev- 
eral facts lead us to use older estimates in form of an- 
alytical approximations: on one hand, Gaidukov experi- 
ments were conducted with a magnetic field intensity of 
H IT, barely around the minimum intensity for the 
FS topology to be relevant, so that a priori we cannot 
expect more than a rough agreement with the experi- 
mental data; on the other hand, dealing with FF known 
only numerically would make so complicated and slow 
our algorithm that it would make sense doing it only 
once it is clear that the algorithm works properly. Since 
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(c) (d) 

FIG. 11: (a) FS of Au according to the Halse formula; (b) nu- 
merical SM of Au in the unit square in coordinates (i/^, Hy, 1) 
obtained using the NTG library; (c) experimental SM for Au 
obtained experimentally by Gaidukov^; (d) numerical SM for 
Au obtained from (b) by symmetry. 

the numerical results proved to be very encouraging, we 
are currently working on adding to our algorithm the 
possibility to work on FF provided in numerical form 
and it is our hope that, now that numerical algorithms 
exist to compare theoretical predictions and experimen- 
tal results, new experiments will be made with stronger 
magnetic fields to allow detailed comparisons. 

Very precise analytical approximations for the FF of 
Au, Ag and Cu have been known since late Fifties; we 
extracted the ones we used from Halse's formulas^^ that 
are known to give rise to surfaces within .2% from the 
actual Fermi Surface: 

£{p)= 3 — ^cos(27rpa^) cos(27rp2y) 
+a(3 - X]cos (47rp^)) 

-\-/3 (3 — XI cos {4:7TPx) COS i^irpy) cos {27TPz)) 

+7 (3 - I] COS {4:7Tpx) COS {4:7rpy)) 

{6 — Y^ (cos {Qirpx) cos (2Tipy) — cos (^npx) cos {Qiipy))) 

a /3 J S Ep 

Cu .00693 -.42501 -.01679 -.03772 1.69167 
Ag -.12030 -.90187 -.14086 -.09483 -0.89789 
Au -.16635 -1.25516 -.09914 -.12704 -2.26213 

Since FS have genus four, in the first Brillouin zone there 
can be at most three cylinders of closed orbits and there- 
fore no more than a pair of "warped planes" , like in the 
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(c) (d) 

FIG. 12: (a) FS of Ag according to the Halse formula; (b) nu- 
merical SM of Ag in the unit square in coordinates {Hx, Hy, 1) 
obtained using the NTC library; (c) experimental SM for Ag 
obtained experimentally by Gaidukov and Alekseevskii^ ; (d) 
numerical SM for Ag obtained from (b) by symmetry. 

genus-3 case, but this time there are six (rather than 
four) basic cycles. Exactly like in the "prison bars" case, 
the symmetries of the FS are such to allow to restrict the 
numerical exploration to half of the square [0, 1] x [0, 1] 
representing the directions H = {Hx^ Hy^l). We per- 
formed numerical simulations for Gold and Silver sam- 
pling the unit square with a lattice of 100 x 100 points. 
The resolution we used for generating a simplicial de- 
composition of the Fermi Surface was again a 100 x 100 
lattice in the cube. The symmetry of the picture with 
respect to the diagonal of the first quadrant is a useful 
extra check about the correctness of the output we get. 

The biggest four zones have been found in both Gold 
and Silver and are labeled exactly by the direction 
at their centers, namely (0,0,1), (1,0,1), (0,1,1) and 
(1, 1, 1) (as it had to be according to Dynnikov^^). 

A fifth zone labeled by (3, 3, 4) has been found in Gold 
only; it does not contain the direction with its label and 
is the only one that has not been detected in Gaidukov's 



experiment. Since this zone is very small it is not clear 
whether it has not been seen because of the limits of 
the semiclassical approximation or of the low intensity 
of the magnetic field used in Sixties by Gaidukov. New 
experimental data taken with a stronger magnetic field 
should be able to answer this question. In fig tTTT d) and 
IT^ d) we extend by symmetry those data and show the 
resulting stereographic projection comparing it with the 
ones generated by Gaidukov from its experimental data. 
In both cases the pictures are in very good agreement. 



VII. CONCLUSIONS 

The method described in this paper allows, for the first 
time, to predict the Stereographic Map of a metal from 
the knowledge of its Fermi Surface. 

Applying this method we first verified, with a toy 
model, Dynnikov's discovery of the fractal structures of 
"global SM" and then produced SM for Au and Ag based 
on analytical expressions for their FS available in liter- 
ature. The zones found numerically for the toy model 
agree perfectly with the corresponding boundaries found 
analytically; the ones found for Au and Ag are very close 
to the corresponding experimental data, that is even 
more remarkable considering that those data, obtained 
more than forty years ago, were taken just at the lower 
threshold for the semiclassical approximation to hold. 

A new zone, not detected experimentally, was found in 
Au. It is not clear whether this is only due to unaccu- 
racies in the FS utlized in calculations or rather to the 
too weak magnetic field used in the Gaidukov experiment 
or possibly to purely quantum mechanical corrections to 
the semiclassical approximation. It is our hope new ac- 
curate Stereographic Maps for Au will be produced to 
help finding out the right answer. 
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